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Abstract. Detailed evolutionary models of the visual bi- 
nary a Centauri, including pre main-sequence evolution, 
have been performed using the masses recently deter- 
mined by Pourbaix et al. (1999). Models have been con- 



structed using the CEFF equation of state, OPAL opac- 
ities, NACRE thermonuclear reaction rates and micro- 
scopic diffusion. A ^-minimization is performed to de- 
rive the most reliable set of modeling parameters p = 
{ta Cen ; «A, cub}, where t a ccn is the age of the sys- 

tem, Y\ the initial helium content, \^]\ the initial metal- 
licity and, a a and ae the convection parameters of the 



two components. Using the basic Bohm-Vitense (1958) 
mixing-length theory of convection, we derive pev = 
{2710 Myr, 0.284, 0.257, 1.53, 1.57}. We obtain a notice- 
ably smaller age than est imate d previously, in agree- 
ment with Pourbaix et al. ( |l999 ), mainly because of the 
larger masses. If convective core overshoot is considered we 
get p ov = {3530 Myr 0.279 0.26 4, 1.64, 1.66}. The use of 
Canuto & Mazitelli ( [1991 , 1992; ) convection theory leads 
to the set p CM = {4086 Myr, 0.271, 0.264, 0.964, 0.986}. 
Using the observational constraints adopted by Guenther 
& Demarque (2000), and the basic mixing- length theory, 
we obtain p GD = {5640 Myr, 0.300, 0.296, 1.86, 1.97} and 
surface lithium depletions close to their observed values. 

A seismological analysis of our calibrated models has 
been performed. The determination of large and small 
spacings between the frequencies of acoustic oscillations 
from seismic observations would help to discriminate be- 
tween the models of a Cen computed with different masses 
and to confirm or rules out the new determination of 
masses. 
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1. Introduction 



As the Sun's nearest stellar neighbors, the two members 
of the visual binary a Centauri A & B (G2V + K1V) 
provide the most accurate potentiality of testing stellar 
physics in conditions slightly different from the solar ones 
and then deserve undivided attention for internal struc- 
ture modeling and oscillation frequencies calculations. By 
coincidence, the masses and the spectral types of compo- 
nents A/B (HD 128620/1, IDS 14328-6025 A/B, Hippar- 
cos 71 683/1) bracket those of the Sun. The high apparent 
brightness and the large parallax imply that surface abun- 
dances and astrometric parameters are known better than 
for any star (except the Sun. The basic intent of this paper 
is to model a Cen A & B using updated physics. This al- 
lows to predict of the p-mode oscillation frequencies, which 
will be useful to exploit future asteroseismological obser- 
vations as expected, for instance, from the project Con- 
cordiastro at the South Pole (Fossat et al. [2000 ). Also, as 
"solar like stars" , the two com ponent s are primary targets 
for the MONS (Kjeldsen et al. |l999bj ) spatial mission and 
their oscillations are expected to be well separated in the 
frequency spectrum. 

Based on the reasonable hypothesis of a common origin 
for both components, i.e. same initial chemical composi- 
tion and age, the calibration of a binary system consists 
in determining a consistent evolutionary history for the 
double star, given (1) the positions of the two compo- 
nents in a H-R diagram, (2) the stellar masses and, (3) 
the present day surface chemical composition. The goal is 
to compute evolutionary models that reproduce the ob- 
servations. This procedure yields estimates for the age i*, 
the initial helium mass fraction Y{ and initial metallic- 
ity [^]i (logarithm of the number abundances of iron to 
hydrogen relative to the solar value), which are funda- 
mental quantities for our understanding of the galactic 
chemical evolution. We also derive values of the "mixing- 
length parameter" or "convection parameter" a, ratio of 
the mixing- length to the pressure scale height. Once the 
initial masses and the physics are fixed, the modeling of 
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the two components A & B of a binary system requires a 
set p of five so-called modeling parameters: 

f Fe 

P = S t*,Y[, [— ]i,aA,a B 

In most cases there are only four reliable observables 
namely, the effective temperatures T c g a, T c ffB and the 
luminosities La, Lb or the gravity log<?A, log<?B of each 
component. Therefore, one of the unknowns has to be 
fixed. Very often the mixing-length parameters are as- 
sumed to be the same for both components, even if the 
mass ratio differs significantly from unity. Once detailed 
spectroscopic analyses have been performed on the sys- 
tem, the precise present day surface metallicities of stars 
come as additional observational constraints. 

In this paper we attempt to reproduce the observed 
metallicities and, if possible, the lithium depletion by 
means of models including microscopic diffusion. We cali- 
brate the binary system using both Bohm-Vitense's ( 



1958 



1992 



hereafter MLT B v) and Canuto & Mazitelli ( |1991 
hereafter MLTcm) mixing- length convection theories. 

The paper is divided as follows: in Sect. ||, we recall 
the main results obtained in previous theoretical works 
and, in Sect. |^, we emphasize the difficulties related to the 
choice of mixing-length parameters. In Sect. we discuss 
the modeling of the transport processes acting beneath 
the convection zone. The observational material, relevant 
to the evolutionary status of a Cen and available in the 
literature, is collected in Sect. ||. The method of calibration 
is described in Sect. [| In Sect. 0, we present the stellar 
modeling procedure. In Sect. |[ we give the results with 
emphasis on the seismological analysis. We summarize our 
results and conclude in Sect. ||. 

2. Previous theoretical works 

On the basis of rotational velocity measurements, ob- 
served Li abund ances and Call emission intensity, Boes- 
gaard & Hagen (1974) attributed to the system an age of 
t a Cen = 3.6 Gyr. Then, several groups calculated stellar 
evolution models to calibrate the system and draw infor- 
mation on the two components as well as on the physics 
governing their structure. Table [j] gives the values of the 
calibration parameters derived in all those studies which 
we now briefly summarize. 

Initially, Flannery & Ayres (1978) and Demarque et 
al. Q1986D could only use the luminosities of a Cen A & 



B as observational constraints. Flannery & Ayres models 
support the fact that the system is metal rich with respect 
to the Sun with Z a Q cn ~ 2Z Q . Demarque et al. (1986) 
derived the age of the system as a function of metallic- 
ity; they also computed the p-mode oscillation spectrum 



of a Cen A. Noels et al. (1991) introduced a general proce 



dure for fitting models to the binary system. They derived 
the age, the helium content and the metallicity of the sys- 
tem and the value of the MLTbv parameter assuming that 



it is the same for the two stars. Neuforge ( 1993| ) revis- 
ited that work using OPAL opacities (Iglesias et al. 1992) 
complemented by her own low-temperature opacities. Fer- 
nandes & Neuforge ( [1995 ) showed that the mixing- length 
parameter a obtained through calibration is different for 
the two stars and that their values become very similar 
if the mass fraction of heavy elements increases above 
Z ~ 0.035. They also performed model calibrations with 
the MLTcm convection treatment with a mixing-length 
equal to the distance to the top of the convective enve- 
lope, thus avoiding the calibration of a c onvection param- 
eter. In parallel, Edmonds et al. ( 1992| ) were the first to 
add the observed metallicity as a constraint (releasing in 
turn the hypothesis that a is unique) and to include the 
effec t of m icroscopic diffusion. More recently, Pourbaix et 
al. ( 1999 ) revisited the calibration of the a Cen system. 
They calculated a new visual orbit on the basis of avail- 
able separations, position angles and precise radial veloci- 
ties measurements and derive new consistent values of the 
orbital parallax, sum of masses, mass ratio and individual 
masses. The main result is that the masses of the com- 
ponents are 5% higher and that the helium abundance 
and age are significantly smaller than previous estimates. 
Guenther & Demarque ( |2000 ) performed several calibra- 
tions of the system using different values of the parallax 
including the Hipparcos value and models calculated with 
updated physics including helium and heavy elements dif- 
fusion. They estimated the uncertainties on the calibrated 
parameters resulting from the error bars on mass, lumi- 
nosity, effective temperature and chemical composition. 

The differences between these calibrations reflect the 
great improvements in the description of the stellar micro- 
physics achieved during the last two decades and the 
progress of the analysis of the observational data. How- 
ever, uncertainties remain on the transport processes (e.g. 
convection, microscopic and turbulent diffusion). Appar- 
ently none of the previous calibrations has attempted to 
reproduce the obs erved surface metallicities (except Guen- 
ther & Demarque 200C ) and lithium abundances which are 
fundamental data for the understanding of the transport 
processes beneath the convection zone of solar-like stars. 

As asteroseismology will in a near future strongly con- 
strain the stellar models, some of the above described the- 
oretical works give the main characteristics of the oscilla- 
tion spectrum of the components (Ai>q and 5vq2 ~ see 
definitions Sect. 8.3.1 ). The "mean" large separation Ai>q 
between the frequencies of modes of a given degree and 
of consecutive radial order depends on the stellar radius 
and mass. The "mean" small separation (5^02 between the 
frequencies of modes with degree I — and 2 and consec- 
utive radial order, measured by (5^02, is sensitive to the 
structure of the stellar core. Table [l] gives estimates of 
these quantities from the literature and from this work. 
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Table 1. Modeling parameters and main characteristics of the oscillation spectrum of aCenA & B derived in this 
study and taken from the literature. The symbols have their usual meaning (see text). The references are: (1) this 
paper, models Abv & Bbv> (2) this paper, models A ov & B ov , (3) this paper, models Acm & Bcm, (4) this paper, 
mo dels A GD & B GD , (5) Flannery & Ayres ( |1978|), (6) Demarque et al. fll986|), (7 ) Noels et al. fll99ip , (8) Edm onds c t 
al. ( |1992| ), (9) Neuforge ( |1993| ), (10) Lydon et al. (l993[ ), (11) Pourbaix et al. ( |1999[ ) , (12) Guenther & Demarque ( |2000D . 
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3. The puzzle of the mixing-length parameter 

The calibration of a solar model provides the value of 
the solar mixing-length parameter a© which in turn is 
currently used to model other stars with the same in- 
put physics. The value of a© changed in time because 
of successive updates of the input physics, in particular 
the low-temperatures opacities and model atmospheres. 
Apart from the difficulty of evaluating a©, the question 
arises whether stars of different masses, initial chemical 
composition and evolutionary status can be modeled with 
a unique a. Because models are also used to calculate 
isochrones, it is important to try to clarify the situation: 
if a is proved to vary significantly with mass, metallicity 
and with evolution, the isochrones might have quite dif- 
ferent shapes which could modify the age estimates. On 
the other hand, models of various masses, all computed 
with a©, can reproduce quite well the slope of the m ain- 
sequence of the Hyades clust er (Pe rryman et al. 1998 ) and 
of field stars (Lebreton et al. 1999 ) observed by Hipparcos, 
indicating that a does not vary much for masses close to 
the solar mass. 

For binaries with well-known properties as a Cen, the 
question of the universality of the mixing-length param- 
eter has been examined many times. Table [l] shows that 
among the attempts of calibrating a Cen A & B in lumi- 
nosities and radii using models based on the MLTbv, two 



yield values of a different for each component and dif- 
ferent from a© (Lydon et al. 1993| , Pourbaix et al. 1999 ) 
while other calibrations suggest similar values for the two 
stars (Edmonds et al. 1992|, F ernandes & Neuforge 1995 ). 
Also, Fernandes et al. ( 1998 ) calibrated three other bi- 
nary systems and the Sun with the same input physics and 
concluded that a is almost constant for [^] in the range 
[^]© ± 0.3 d ex and masses in the range 0.6 — 1.3 M Q while 
Morel et al. ( 20001: ) found a small difference of a (« 0.2) 
in the two components of the l Peg system. 

As pointed out by many authors (e.g. Lydon et 
al. 1993, Andersen 1991) accurate masses and effective 
temperatures are required to draw firm conclusions on 
whether the MLTbv parameters are different or not. This 
is also illustrated by the large error bars on a (w 0.8) 
quoted by Pourbaix et al. (1999) which include all possi- 
ble observable sources of errors including the errors in the 
masses. For a detailed estimate of the error budget see e.g. 
Guenther & Demarque ( |2000 ). 

Some 2-D and 3-D numerical simulations of convection 
have been performed and translated into effective mixing 



length parameter (Abbett et al. 1997, Ludwig et al. 1999, 
Freytag et al. |1999 ). They suggest that, for stars with 



effective temperature, gravity and metallicity close to so- 
lar the mixing-length parameter remains almos t con stant. 
Similar results are obtained by Lydon et al. ( 1993 ) who 
compared models of a Cen A & B based on the results of 
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numerical simulations of convection and found that a sin- 
gle value of a may be applicable to main sequence solar 
type stars. In the MLTcm convection theory, Canuto & 



Mazitelli (1991, 1992) used two modern theories of turbu- 
lence to establish a formula for the turbulent convective 
flux which replaces the MLTbv expression (because the 
full spectrum of convective eddies is taken into account, 
the convective efficiency is magnified by about one order of 
magnitude) . The MLTcm theory also makes use of a mix- 
ing length which is equal, either to the distance towards 
the outermost limit of the convection zone, or to the frac- 
tion (a) of the pressure scale height. The solar calibration 
with the MLTcm treatment yields an a Q value of the or- 
der of unity. It is to be noted that for the Sun, models 
calibrated with MLTcm have frequencies which are closer 
to the observations than tho se ca librated with MLTbv 
(Christensen-Dalsgaard et al. 



1996|) 



4. Diffusion and transport processes 

For the Sun, it is now well established that microscopic 
diffusion by gravitational settling has a measurable im- 
pact on models thus on their pulsation spectrum. As time 
goes on, all chemical species sink in the gravitational well, 
except the lightest one, 1 H. At the surface, the amount 
of hydrogen is enhanced while helium and heavy elements 
are depleted. This implies a decrease of the surface metal- 
licity [Ep] s with respect to time. Since the time-scale of 
abundance variations in the convection zone is roughly 
proportional to the squ are ro ot of the mass of the con- 
vection zone (Michaud 1977 ), the microscopic diffusion 
time-scale decreases with increasing stellar mass imply- 
ing that diffusion is more efficient in massive stars. As a 
consequence, a smaller surface metallicity is expected for 
a Cen A than for a CenB. This is an additional constraint 
for the models. 

As lithium burns at temperatures close to 3MK, the 
observed surface abundance [Li] s (H = 12) is often used 
as a probe for theories of transport processes at work be- 
neath the external convective zone. The microscopic dif- 
fusion alone is not efficient enough to explain the lithium 
depletion observed in the Sun and stars; thus, there is 
a need for an unknown physical process, e.g. a turbulent 
mixing, acting in the radiative zone just beneath the outer 
convection zone - e.g. Schatzman (1996), Montalban & 
Schatzman fl2000p for a review. According to the present 
state of the art, either the shear resulting from the dif- 
ferent rotational status between the convection zone and 
the radiative interior (e.g. Zahn 1992| ), or internal waves 
(e.g. Montalban & Schatzman 200C| ), are believed to be 
responsible for that turbulence. For stars more massive 
than the Sun, turbulent mixing is required, in some (hy- 
pothetical) mixed stellar mass fraction, to connect and 
to extend somehow the mixing of the external convec- 
tion zones. This avoids the complete segregation of helium 
and heavy elements from the surface and can explain the 



AmFm phenomenon (e.g. Richer et al. 2000). Despite nu- 
merous attempts, there is today no fully satisfactory pre- 
scription for turbulent diffusion able to account for the 
observed solar lithium depletion; for other stars the sit- 
uation is even worse ! Further information on the mixing 
processes at work below the outer convection zones can 
be provided by simultaneous observations of lithium and 
beryllium, this latter being depleted deeper in the stars at 
temperatures slightly larger than for lithium. 

5. Observations of the visual binary a Cen 

The calibration of a binary system is based on the ad- 
justment of the stellar modeling parameters, so that the 
model of each component at time t a c on reproduces the 
available observables within their error bars. Among dif- 
ferent possibilities, one has to choose the most suited set 
of variables to be used in the calibration according to the 
observational and theoretical material available. As afore- 
mentioned, because of the assumption of a common origin 
of both components, five unknowns enter the modeling 
of their evolution. There are six most currently available 
constraints: (1) the effective temperatures for a Cen A & 
B, derived from precise detailed spectroscopic analysis, (2) 
either the spectroscopic gravities or the luminosities, these 
latter derived from the photometric data and parallax, us- 
ing bolometric corrections (3) the surface metallicities. In 
the following, we discuss all the relevant available obser- 
vations and we choose the most appropriate and accurate 
observables to be fitted in the calibration process. 

5.1. Astrometric data 

For a visual binary, the parallax uncertainty is the 
largest source of error in the derivation of the sum of 
masses via the third Kepler law (Couteau 1978, 40, 
VI). Unfortunately, a Cen was "difficult" for Hipparcos 
with the secondary at the edge of the sensitivity pro- 
file (Soderhjelm 200C). Soderhjelm improved the adjust- 
ment of the Hipparcos parallax by taking into account the 
known orbital motion but, as expected, the orbital curva- 
ture covered during the 3.25 yr mission, was not sufficient 
to derive a reliable mass ratio of the components. Table 
shows that the new parallax value differs from the former 
by more than 2a. With respect to the original Hipparcos 
value this new one is closer to the long focus photographic 
determin ation of Kamper & Wesselink (1975) based on 
Heintz's ( 1958 ) orbital elements. 



After numerous independent investigations during al- 
most 100 years, Pourbaix et al. ( 1999] ) derived an orbit 
based both on visual astrometric and precise radial veloc- 
ity data, leading to consistent determinations of orbital 
parameters, orbital parallax, sum of masses, mass ratio 
and, thus, precise mass values for a Cen A & B. Here, 
we prefer to use of P ourba ix et al. ( 199S ) mass values 
because Soderhjelm's ( 200fj|) ones are based on disparate 
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Table 2. Astrometric properties of the aCentauri A & B binary system. P is the orbital period in yr, a the semi 
major axis in arc second, w the parallax in mas, K = Mb/(Ma + Mb) the fractional mass. 



K 



M a /Mq M B /M e 



References 



79.92 



17.515 



79.90 ±0.01 17.59 ±0.03 



750 

749 ±5 
750.6 ±4.6 
742.1 ± 1.4 

737.0 ± 2.6 

747.1 ± 1.2 



0.453 

0.454 ± 0.02 



0.45 ±0.1 



1.09 
1.10 



0.90 
0.91 



1.16 ±0.03 0.97 ±0.03 



Heintz ( [19581 ), ( |l982j ) 
Kamper & Wesse link (|l978| ) 
Dema rque et al. ( 1986 ) 
ESA fll997| ) 



Pourbaix et al. (h 999| ) (orbital parallax) 
Soderhjelm (]2000| 



sources (i.e., Hipparcos parallax, Hcintz's (1958) orbit and 
the mass ratio of Kamper & Wesselink ( 1978| )). We do not 
take into account mass errors in the calibrations. 

The inclinations of th e orbits of Heintz (|l958t |l982|) 
and Pourbaix et al. ( 1999 ) agree within less than one tenth 
of a degree, i = 79°29 ± 0°01, sini ~ 1 (the Sun lies close 
to the orbital plane of aCen). With the reasonable as- 
sumption that the orbital and rotational axis are parallel, 
owing to the high inclination, a Cen A & B are seen about 
equator on and the observed rotational velocities are close 
to their equatorial values, i>sini ~ v. 



5.2. Photometric data 

Photometric data are available in the V, B (Mermilliod 
et al. |1997D and K (Thomas e t al. 



Table 

& B. 



1973) bands. But, as 



claimed by Flannery & Ayres (1978) and Chmielewski et 
al. ( |l992p , the dynamical range of commonly used pho- 
tometers is limited, which does not allow precise pho- 
tometric measurements of bright stars. This makes the 
a Cen photometric data rather unsafe (Chmielewski et 
al. |j"992] ) and we preferred to use the spectroscopic data. 
However, for comparison, we have used these photometric 
data to derive bolometric corrections and effective tem- 
peratures of the two compo nents from (1) the empirical 
calibrations of Alonso et al. ( 1995| , 1996) a nd (2) the the- 
oretical calibrations of Lejeune ct al. ( 1998 ). T he T c k (V — 
K, [^]) empirical calibration of Alonso et al. (|l996| ) leads 
to Toff, A = 5738 ± 40 K andT Gf f iB = 5165 ± 70 K while the 
theoretical (ba sed on model atmospheres) calibrations of 
Lejeune et al. (|l998|) yield T cS . A ~ 5800 K and T cS . B ^ 
5250 K, all well compatible with spectroscopic tempera- 
tures. The Fb \(V — K, [^]) e mpiri cal calibration of the 
bolometric flux (Alonso et al. 1995 ) provides bolometric 
corrections BC A {V) = -0.09 and BC B (V") = -0.23 close 
to the values i nferre d from the theoretical calibrations of 
Lejeune et al. (199?) (—0.11 and —0.23 respectively). Let 
us point out that these bolometric corrections are larger 
than those adopted by Guenther & Demarque (200C). 



4. 

The 
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Thevenin ( |1998D 



Atmospheric abundances of a Cen A 
references are: (1) Furenlid & Mey- 
Neuforge-Verheecke & Magain ([1997D, (3) 
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0.05(1) 




0.14 (2) 


Zr 


0.20 (3) 


0.17(2) 


0.04(1) 






La 


0.40 (3) 










Ce 


0.40 (3) 










Nd 


0.35 (3) 


0.20 (2) 








Sm 


0.40 (3) 










Eu 




0.15(2) 









5.3. Spectroscopic data and atmospheric parameters 

For both stars many spectroscopic data exist in the lit- 
erature with good signal to noise ratios because of their 
high magnitudes. As Table || exhibits rather scattered re- 
sults we decided to reestimate the effective temperatures, 
surface gravities and metallicities. For our analysis we 
used the spectroscopic data published by Chmielewski ct 



al. (1992) and Neuforge-Verheec ke fe Magain (1997), the 



atmosphere models of Bell et al. ( 1976 ) grid and the basic 
technique of the curve of growth. 



6 



P. Morel, J. Provost, Y. Lebreton, F. Thevenin and G. Berthomieu: Calibrations of a Cen A & B 



Table 3. Spectroscopic data of the a Cen system. The effective temperatures are in K, the metallicities and the 
lithium abundance in dex. For sake of brevity the uncertainties are not recalled. 



TeffB loggA loggB [] 



]b [Li] a [Li] 



References 



5793 


5362 






0.22 


0.12 


5727 


5250 


4.38 


4.73 


0.29 


0.37 


5820 


5350 






-0.01 


-0.05 


5720 




4.27 




0.15 




5793 


5305 


4.40 


4.65 


0.20 


0.20 


5793 


5305 


4.50 


4.50 


0.20 


0.14 


5727 


5250 


4.42 


4.65 


0.20 


0.26 


5727 




4.27 




0.10 




5800 


5325 


4.31 


4.58 


0.22 


0.26 


5720 




4.27 




0.15 




5830 


5255 


4.34 


4.51 


0.25 


0.24 


5730 


5250 


4.2 


4.6 


0.20 


0.22 



1.28 



1.3 



1.37 



< 0.4 



French & Powe ll ( p.971p 
England (|l973j ) 
Bessell ( |L98lD 

Soberblom & Dravins (|l984j) 
Rabolo et al. (|l986[ 
Smith et al. ( | l98 
Abia et al. (Il98i 




Edvardsson Ql988| ) 
Furenlid & Meylan (|1990| ) 
Chmielewski et al. (|1992|) 

Edvardsson et al. (11993) 

Neuforge-V erheec ke & Magain (1997) 
King et al. (|l997|) 
Thevenin (|1998j) 



5790 5260 4.32 4.51 0.20 0.23 



This paper 



5.3.1. Effective temperatures 

Table || presents effective temperatures compiled from the 
literature. Chmielewski et al. ( 1992 ) derived the effective 
temperatures from the best fit of the wings of the H Q line. 
Recently Neuforge-Verheecke & Magain (1997) obtained 
the T c ff by forcing the abundances derived from Fel lines 
to be independent of their excitation potential. They not 
only used the Fel lines to determine the effective temper- 
ature of each star but also checked their results using the 
H Q line profiles. In a Cen A, the two methods lead to con- 
sistent results that are also i n good agreement with those 
of Chmielewski et al. (1992), although the two analyses 
made use of different model atmospheres. In aCenB, the 
two methods lead to marginally discrepant results possi- 
bly due to Non-LTE effects and to the sensitivity of H Q 
to the treatment of convection. 

We have reconsidered the fit of hydrogen lines of 
Chmielewski et al. (1992). We took into account the over- 
abundance of the a Cen system which strengthens the 
metallic lines and displaces the pure wings of H lines; we 
also took into account the relative strength of metallic 
lines in aCenB, different from what is found in solar- 
like profiles, because of the lower temperature. Conse- 
quently, the fit proposed by Chmielewski et al. (1992) 



seems to overestimate the temperature and we displace 
it toward cooler values, more importantly in the case of 
a Cen B. Our adopted effective temperatures are respec- 
tively 5790 ± 30 K and 5260 ± 50 K for a Cen A & B. These 
new proposed temperatures agree, within the 1 a error 
bars, with the Fel based temperature determinations of 



Neuf orge -Verheecke and Magain ( |1997| ) and, as quoted in 
Sect. 5^, they are close to the values based on Lejeune et 
al. ( 1998 ) calibrations. 



5.3.2. Surface gravities 

Usually logg, the logarithm of the surface gravity, is de- 
termined by the detailed spectroscopic analysis using the 
ionization equilibrium of iron. Table [5] shows the results 
obtained by different authors. Wc remind that the surface 



gravities from Neuforge-Verheecke & Magain (1997) and 
from Chmielewski et al. (1992) are very close. Using the 
published equivalent width (hereafter EW) and the error 
bars of Neuforge-Verheecke & Magain, we derive by the 
basic curve of growth procedure log^A = 4.32 ± 0.05 and 
log 3b = 4.51 ± 0.08 for a Cen A & B, respectively. 

5.3.3. Chemical composition 



Furenlid & Meylan (199C) and Neuforge-Verheecke & 
Magain (1997) have found an average metal overabun- 
dance of around 0.2 dex. As in Sect. 5.3.2) , using oscilla- 
tor strengths, microturbulence parameter, EW, published 
by Neuforge-Verheecke & Magain and the iron curve of 
growth derived as aforesaid, we get [§]a = 0.20±0.02dex 
and [§]b = 0.23±0.03dex for aCenA & B, respectively; 
the error bars are from Neuforge-Verheecke & Magain. We 
note that metallicities agree, within the 1 a error bars, 
with the values of Neuforge-Verheecke & Magain. 

Table |] shows that C, N, O and Fe, the most abun- 
dant "metals" and electron donors, have uniform over- 
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abundances of about ~ 0.2 dex. In consequence, it is rea- 
sonable to construct models having a uniform abundance 
of all metals enhanced by a factor 1.7, corresponding to 
0.2 dex compared to the Sun, and use opacities with solar 
mixture. On the main-sequence, for stellar masses close 
to the solar one, the differential segregation between the 
metals does not significantly modify the models (Turcotte 
et al. 1998; ) this permits, for the equation of state and 
opacity data, to keep the ratios between "metals" to their 
initial value despite the differential segregation by micro- 
scopic diffusion and gravitational settling. 

Table || reveals a good agreement between the mea- 
surements of the lithium abundance. aCenA presents a 
depletion close to the solar one ([Li] S Q = 1.11, Grevesse 
& Sauval 1998) while, in a CenB, the lithium is clearly 
more depleted. The assumption of a common origin of the 
two stars eliminates the possibility of a different initial 
lithium abundance in the pre-stellar nebulae. Neither an 
error on temperatures nor on surface gravities in model 
atmospheres can explain these differences with the solar 
depletion. 

The situation is similar for 9 Be. Observations by Pri- 
mas et al. ( 1997| ) indicate that 9 Be is depleted in a CenB 
while aCenA exhibits a quasi-solar abundance. 

5.3.4. Rotation 



found an evidence for almost the same periodicity as that 



For aCenA, Boesgaard & Hagen (|l974j) estimated a ro- 
tational period 10% la rger than the solar one. More re- 
cently Saar & Osten ( 1997 ) measured rotational veloci- 
ties for aCenA & B of respectively 2.7 ± 0.7 km s -1 and 
1.1 ± 0.8 km s -1 , in good agreement with other determi- 
nations. With the estimated radius of i? Q CenA — l-2i? Q 
and fl Q cenB — O.91i?0, the periods of rotation are respec- 
tively P Q cenA — 22 d, and P a c on B — 41 d. They bracket 
the solar value. 

5.4- Seismological observations 

Many attempts have been made to detect solar-like p- 
mode oscillations in a Cen A, from ground-based observa- 
tions. They lead to controversial results. Gelly et al. (1986) 
reported a detection, but the claimed mean large separa- 
tion between the oscillation frequencies Ai>o = 165.5 pHz 
was very large and inconsistent with the theoretical value 
(Demarque et al. 1986). Despite the use of more sensi- 
tive techniques, Brown & Gilliland (1990) failed to detect 
any oscillation and gave an upper limit for the ampli- 
tude of 0.7ms" 1 , i.e. 2-3 times the solar one. Pottasch 



et al. (1992) have detected oscillations in a Cen A. They 
found a set of regularly spaced peaks, with a mean sep- 
aration corresponding to Ai>o — 110 /xHz, consistent with 
the known properties of this star, but with an ampli- 
tude of oscill ation larger than the upper limit of Brown 
& Gilliland ( |l990| ). Edmonds & Cram (|l995j ) did not 
detect unambiguously the p-mode oscillations, but they 



found by Pottasch et al. (1992) in the power spectrum. 
More recently, Kjeldsen et al. (1999a) found tentative ev- 
idence for p-modc oscillations in aCenA, from equiva- 
lent width measurements of the Balmer hydrogen lines. 
They proposed four different possible identifications of the 
eight observed frequency peaks in the power spectrum, 
which correspond to sets of the large and small frequency 
spacings Av and ^ ,2 of (106.94,12.30), (106.99,8.15), 
(100.77,11.70) and (100.77,6.42), in ^Hz. The frequency 
resolution of the seismological observations will be much 
improved in a near future, up to 0.1 /zHz, with several 
ground-based and space projects. 

5.5. Adopted observables constraining the models 

To constrain the models at present day, contrarily to pre- 
vious works, we prefer to use the spectroscopic gravities 
we derived consistently from effective temperature and 
metallicities instead of the luminosity derived from the 
photometry, bolometric correction and parallax. Table 
gives the effective temperatures, gravities and metallicities 
we selected to constrain the models. The masses are from 



Pourbaix et al. (199G). We emphasize that the parallax 
is required to derive the masses but neither the effective 
temperatures nor the gravities. 

While t his wo rk was under investigation, Guenther & 
Demarque ( 200C ) published new calibrations of the a Cen 
binary system. As a matter of comparison we have cal- 
ibrated the a Cen binary system with their constraints. 
Table || presents the observable con strain ts of the "best" 
model of Guenther & Demarque ( 120001) . They signifi- 
cantly differ from ours, though correspondingly to al- 
most the same mass ratios. The trigonometrical paral- 
lax of Sodcrhjelm ( 200C) is larger than the orbital paral- 
lax of Pourbaix et al. ( 1999 ) and leads to smaller masses 
and larger (respt. smaller) luminosity for a Cen A (respt. 
a CenB). 



6. The method of calibration 

The calibration of a binary system is based on the ad- 
justment of stellar modeling parameters to observational 
data at the age of the system. For a given mass, fixing 
the physics, the effective temperature, surface gravity and 
metallicity of a model have the formal dependences with 
respect to modeling parameters: 



Jeff (*)mod 
logffWmod 

Fe 

["jr"]s(*)mod 



Fe 



Fe / Fe 
[— ] s lU;Y u [— ] i; a* 
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The basic idea of the x 2 fitting has been developed by 
Lastennet et al. (1999). To find a set of modeling param- 
eters: 



Pa Cc 



v Fe. 
iaCcn, ii, [ jj Ji, «A, C*B 



leading to observables as close as possible to the obser- 
vations T c ff(*), logg(*) and [§] s (*), we minimize the 
X 2 (t*, Yi, ctA, cub) functional defined as: 

2 



X 



E 

*=A,B 



Tcff (*)mod — ^cff(*) 

<r(T off (*)) 



+ 



+ 



(1) 



/ logg(*)mod - logg(*r 
V o-(logsW) 

/ [f] B (*)mod-[f]s(*) \ 2 
l *([§].(*)) j ' 

where cr (T c ff (*)), cr (log <?(*)) and er ([§] s (*)) are the un- 
certainties associated to each star. For a grid of modeling 
parameters p, we have computed main-sequence evolution 
of models with a Cen A & B masses. Then the x 2 was com- 
puted using Eq. (jlj) in a refined grid obtained by interpo- 
lations. We kept for the solution the "best" p = p a Cen 
which corresponds to the Xmin- These best modeling pa- 
rameters are used to compute models of a Ccn A & B 
including pre main-sequence and their frequencies. We 
do not further attempt neither to improve nor to inves- 
tigate the stability of the solut ion via the steepes t de- 
scent method (Noels et al. |l99lj Morel et al. |2000b|) . Ta- 
ble [j] shows the confidence limits of modeling parameters 
of models computed in this paper. The confidence limits 
of each modeling parameter, the other being fixed, corre- 
spond to the maximum/minimum values it can reach, in 
order that the generated models fit the observable targets 
within their error bars. 



7. Evolutionary models 

Basically the physics employed for the calculation of mod- 



els is the same as in Morel et al. (1997). The ordinary as- 
sumptions of stellar modeling are made, i.e. spherical sym- 
metry, no rotation, no magnetic field, no mass loss. It has 
been already claimed, see Morel et al. (2000a), that stellar 
models of about one solar mass, computed with the same 
physics but initialized either on pre main-sequence or at 
zero-age homogeneous main-sequence are almost identical 
after a few 10 7 years, a small quantity with respect to the 
expected age of the a Cen system. To save computations, 
along the search of a solution with the \ 2 minimization, 
we have initialized each evolution with a homogeneous 
zero-age main-sequence model. The models presented in 
Table |B| include the pre main-sequence evolution. They 
are initialized with homogeneous zero-age stellar model in 
quasi-static contraction with a central temperature close 
to the onset of the deuteron burning, i.e. T c ~ 0.5 MK. 



7.1. Nuclear and diffusion network. 

We have taken into account the most important nuclear re- 
actions of PP+CNO cycles (Clayton, |1968|) . The relevant 
nuclear reaction rates are taken from the NACRE compila- 



tion (Angulo et al. 1999) with the reaction 7 Be(e 



'Li 



taken from the compilation of Adelberger et al. (1998). 
Weak screening (Salpeter 1954) is assumed. We have used 
the meteoritic value (Grevesse & Sauval 1998) for the ini- 
tial lithium abundance, [Li]; = 3.31 ± 0.04. For the cal- 
culations of the depletion, lithium is assumed to be in its 
most abundant isotope 7 Li form. The initial abundance of 
each isotope is derived from isotopic fractions and initial 
values of Y = 3 He + 4 He and Z in order to fulfill the basic 
relationship X + Y + Z = 1 with X = X H + 2 H. For the 
models computed for the % 2 fitting, the isotopes 2 H, 7 Li, 
7 Be are set at equilibrium. 

Microscopic diffusion is described by the simplified for- 



malism of Michaud & Proffitt (1993) with heavy elements 
as trace elements. We have neglected the radiative acceler- 
ations as they amount only to a tiny fraction of gravity in 
the radiative part f or sta rs with masses close to the solar 
one (Turcotte et al. 1998| ). We assume that changes of Z, 
as a whole, describe the changes of metals and we use the 
approximation: 



Fe 
[— ] 



i i Zfc • 
log(^, 

log(f) ■ 



log( 



Z_ 

X 



i i Zfc \ 
lo g(^)© 



!og(j)0 



with -Hp is the iron mass fraction within Z. 



7.2. Equation of state, opacities, convection and 
atmosphere. 

We have used the CE FF eq uation of state (Christensen- 
Dalsgaard fe Da ppen 1992 ) and the opacities of Iglesias 
& Rogers ( 1996 ) compl emen ted at low temperatures by 
Alexander & Ferguson ( 1994 ) opacities for the solar mix- 
ture of Grevesse & Noels ( 1993 ). We have not taken into 
account the changes of abundance ratios between the met- 
als within Z due to diffusion; for stellar masses close to the 
solar one they do not really affect the structure of models 
(Turcotte et al. |199SP . 

In the convection zones the temperature gradient is 
computed according to either MLTbv or MLTcm convec- 
tion theories. The mixing- length is defined as I = aH p , 
where H p is the pressure scale height. The convection 
zones are mixed via a strong full mixing turbulent dif- 
fusion coefficient d m — 10 13 cm 2 s _1 which produces a ho- 
mogeneous composition (Morel |l997 ). 

At the end of the pre main-sequence both components, 
have for a few million years, a temporary convective core. 
For a Cen A, slightly more massive than the Sun, a second 
convective core is formed during the main-sequence due 
to the onset of the CNO burning (e.g. Guenther & De- 
marque 2000). Following the prescriptions of Schaller et 
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al. (1992) we have calibrated models with overshooting of 
convective cores over the distance O v = 0.2 mm(H p , R co ) 
where R co is the core radius. 

The atmosphere is restored using a grid of T(t, T e g, g) 
laws, provided by Cayrel ( [2000 ), (r is the Rosseland op- 
tical depth) derived from atmosphere models with the 
solar mixture of Grevesse & Noels ( 1993| ) and metallic- 
ity [§] = 0.2. The atmosphere models were computed 



with the Kurucz ( |199l[ ) ATLAS12 package. The connec- 
tion with the envelope is made at the optical depth i\, = 20 
where the diffusion approximation for radiative transfer 
becomes valid (Morel et al. 1994). A smooth connection 



of the gradients is insured between the uppermost layers 
of the envelope and the optically thick convective bot- 
tom of the atmosphere. It is an important requirement 
for the calculation of eigenmode frequencies. The radius 
i?* of any model is taken at the optical depth where 
T(t*) = T e ff. Typically, t* increases from ~ 0.43 in 
the initial pre main-sequence model, until t* ~ 0.53 at 
the present time. The mass of the star M+ is defined as 
the mass enclosed in the sphere of radius R+. The exter- 
nal boundary is located at the optical depth r oxt = 10~ 4 , 
where the density is fixed to its value in the atmosphere 
model /o(r ex t) = 2.22 10 _9 gcm -3 . To simplify, the chemi- 
cal composition within the atmosphere models is assumed 
to be unaffected by the diffusion. 



7.3. Numerics. 

Models have been computed using the CESAM code 
(Morel 1997). The numerical schemes are fully implicit 
and their accuracy is of the first order for the time and 
third order for the space. Each model is described by about 
600 mass shells, this number increases up to 2100 for the 
models used in seismological analysis. Evolutions are de- 
scribed by about 80 models. About half of them concerns 
the pre main-sequence evolution. 



8. Results 

For different sets of modeling parameters we have com- 
puted evolutionar y tra cks using the convection theories 
of Bohm-Vit ense fll958j MLT BV ) and Canuto & Mazitelli 



1991| , |1992j , MLT CM ). Using the x 2 minimization de- 
scribed in Sect. ||, we have selected the best fits taking 
into account the accuracy of the observational constraints. 
Tables [l] and || (bottom panel) present the characteristics 
of a Cen A & B and Sun models computed with the same 
input physics. Figure [l] presents the evolutionary tracks 
in the logL/L© versus logT e ff H-R diagram. We also per- 
formed x 2 minimization with the constraint a a = ob, but, 
for brevity sake, we do not report the results as they do not 
significantly differ from those obtained without this con- 
straint. For all a Cen A models a convective core, caused 
by an enhancement of the CNO nuclear energy genera- 
tion, appears as soon as the central temperature becomes 



slightly larger than 17 MK and the mass fraction of hy- 
drogen X c < 0.37. Figure [j], Tables |l] and || clearly show 
that various sets of modeling parameters allow to verify 
the observational constraints within the error bars. In ab- 
sence of seismological observations there is no criterion 
available to discriminate between models. Tables [l] and || 
show that, within the confidence domains, the values of 
the mixing length parameter of the two components are 
very close. 

8.1. MLTbv models. 

Models Abv & Bbv have an age i Q cen = 27 10 My r 
compatible with the estimate of Pourbaix et al. ( 1999 ). 
The values of the convection parameters are almost equal 
(a a — Qb) and close to the value derived if they are forced 
to be identical in the x 2 minimization. They are quite 
close to the value obtained by Noels et al. ( 1991 ) b ut dif- 
fer notably from the estimate of Pourbaix et al. ( 1999 ) 
though remaining within their large confidence domains. 
They are smaller than the solar value a Q . In a Cen A a 
convective core is formed at time t ~ 2.6 Gyr a few Myr 
just before t Q Cen- Figure [j], bottom panels, shows the loci 
in the H-R diagram reached by the Abv & Bbv models 
when one modeling parameter changes within its uncer- 
tainty domain, the other modeling parameters being fixed. 
The limited extents of the permitted solutions are conse- 
quences of limits in metallicity. Table |^ shows the partial 
derivatives of the observed parameters with respect to the 
modeling ones. Our results are similar to those of Brown 
et al. (1994), except for the derivatives with respect to 



the age which are very dependent on the exact evolution- 
ary stage on the main s equen ce. As already noticed by 
Guenther & Demarque (200C), the luminosities are very 
sensitive to the helium mass fraction (see the large values 
of their derivatives with respect to Yi in Table ||) , which in 
turn gives a narrow confidence level for the initial helium 
abundance (see Table |l|) . 

For both components the surface lithium depletion is 
not sufficient enough to fit the observed values. 

Models A ov & B ov with overshooting of convective 
cores have an age significantly larger than models Abv 
& Bbv- Indeed their outer convective zones penetrate 
deeper. For model B ov this depletion is marginally com- 
patible with the observation, while it is not large enough 
for model A ov . All these differences with models Abv 
& Bbvj of course, result from the overshooting of con- 
vective cores but in a roundabout way via the % 2 min- 
imization which adjusts the modeling parameters as a 
whole. The convective core of model A ov is formed at time 
t ~ 3.16 Gyr. 

Models Aqd & Bgd are calibrated according to the 
Guenther & Demarque ( 2000| ) observational constraints 
and masses. The age, initial helium mass fraction and 
metallicity and the convection parameters are significantly 
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Fig. 1. Evolutionary tracks in the H-R diagram of models Abv, Bbv (full), Acm, Bcm (dashed) and A ov , B ov (dot- 
dash-dot). Dashed rectangles delimit the uncertainty domains. Top panel: full tracks from PMS. The stellar evolution 
sequences are initialized on the pre main-sequence soon after the deuteron ignition. The "+" denote 1 Gyr time 
intervals along the evolutionary tracks. Middle left and right panels: enlargements around the observed aCenA, & 
B loci. Bottom left and right panels: loci of aCenA & B models computed with the MLTbv convection theory and 
modeling parameters within the confidence domains presented in Table [|; full triangle: i Q con changes only, full star: 
a changes only, empty circle: Yj changes only and full dot: changes only. 
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larger than those derived using the mass values of Pour- 
baix et al. (1999). We have obtained similar modeling 
parameters from the x 2 minimization using models with 



Gucnthcr & Dcmarque (200C) masses and our observa- 
tional constraints. Therefore, the differences between the 
modeling parameters of models Abv & Bbv and Agd & 
Bqd mainly result from the mass differences, but also from 
the disparity of observing targets, although in a less ex- 
tent. For models Agd & Bqd the convection parameters 
are close to the solar value. In the model Agd the con- 
vective core is formed at the time t ~ 3.7 Gyr, larger than 
in other MLTbv models and the lithium depletion at the 
surface is marginally compatible with the observed value. 
The value of the central hydrogen mass fraction is half 
that obtained in the other MLTbv models. In model Bgd 
the lithium depletion is compatible with the observation. 

8.2. MLTcm models. 



In agreement with Canuto & Mazitelli (1991, 1992) con- 
vection theory, the convection parameters of models Acm 
& Bcm are close to unity and close to the solar value. 
The age is slightly smaller than the solar one but larger 
than that of models Abv & Bbv! the outer convection 
zone is deeper. According to the observations, the lithium 
of model Bcm is almost totally depleted but practically 
no depletion occurred in models Acm and the solar one. 
The convective core of model Acm is formed at time 
t ~ 3.9 Gyr. 

8. 3. Seismological analysis of a Cen A & B 

The stars a Cen A & B are solar-like stars. The oscillations 
of such stars may be stochastically excited by the convec- 
tion as in the case of the Sun. The amplitudes of solar- 
like oscillations have been estimated for stars of different 



masses and ages (Houdek |1996| Houdek et al. [1999J) . The 
models we have calibrated are close in the H-R diagram, 
but have different internal structure and they could be 
discriminated with the help of seismology. The properties 
of the stellar oscillations are related to the variation along 
the radius of the sound speed c and of the Brunt- Vaissala 
frequency N (Christensen-Dalsgaard & Berthomieu 1991 ). 
The p-mode frequencies are mainly related to c, while the 
(/-modes, with lower frequencies, are determined essen- 
tially by N . We have computed for all our models the adi- 
abatic frequencies of oscillation for modes of low degrees 
£=0 to 3, which may be detected by future observations. 
Table || presents the p-mode frequencies for the models 
Abv and Bbv- As in the solar case (Provost et al. 2000| ). 
the lowest frequency p-mode oscillations of a Cen A have 
a mixed character between p- and <?-mode behavior and 
are very sensitive to the structure of the central stellar 
regions. 



Table 7. Theoretical global asymptotic characteristics 
of the low degree p-mode and (/-mode spectrum of the 
star a Cen A & B. The quantities vqo, Avq Svq^ and 5v\^ 
(no = 21), describing the p-mode oscillations, are given in 
/iHz - see definitions Sect. |S.3.1| . The characteristic (/-mode 
period Pq is given in mn. The lower panel presents the vari- 
ations A[Az/o],A[<52/o,2], A[<5i/i,3] and A[Po] for the models 
computed with the MLTbv theory and extreme modeling 
parameters a, i Q ccn, Y\ and [w]; within the confidence 
domains presented in Table |l|. 
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Po 




Abv 


2391.0 
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8.3.1. p-mode oscillations 

The frequencies of p-modes of given degree £ are al- 
most quasi-equidistant. Figure |^ represents the variation 
of Ais r ^i = v n j — v n ~i.z as a function of the frequency 
for models Abv and Bbv! n is the radial order of the 
mode. Av rh i represents the large spacing between the 
mode frequencies. It is roughly constant at large frequen- 
cies, larger than 2300 /iHz for a Cen A and 3500 /iHz for 
aCenB. Below these frequencies it varies within 10%. In 
the lower frequency range, Av n ^ has an oscillatory behav- 
ior which is the s ignature of the helium ionization zone 
(e.g. Gough 1991 ). We note that it is important to take 
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Fig. 2. Variations of the large frequency separations be- 
tween modes of consecutive radial order Af„^ = v n j — 
Vn-\,i for p-modes of degree I = (full point) and I = 1 
(open point) as a function of the frequency for aCenA 
& B (models Abv & Bbv)- The relation between the fre- 
quencies and the radial order are taken from Table |l As 
asymptotically predicted, Ai/ n i is almost constant at high 
frequency. 



into account such a variation while searching peak equidis- 
tancy in the p-mode power spectrum. 

Another characteristic quantity of the oscillation spec- 
trum is the small separation, i.e. the difference between 
the frequencies of modes with a degree of same parity and 
with consecutive radial order: 

This small quantity is very sensitive to the structure of the 
core mainly to its hydrogen content. Asymptotic analysis 
predicts that it is proportional to t(l + 1). Figure ^ gives 
5^0,2 arid | x Svx,3 as a function of the frequencies for the 
models Abv and Bbv- For a given model, these two quan- 
tities are very close at high frequencies, as expected from 
asymptotic approximation and they vary almost linearly 
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Fig. 3. Variations of the small frequency differences <5fo,2 
and | <5^i, 3 as a function of the frequency for a Cen A & B 
(models Abv & Bbv)- Same symbols as in Fig. g. 
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Fig. 4. Variations of the differences of frequencies A n ^ = 
v n £ — 2v n+ ij + v n +2,e for p-modes of degree £ = 0,1,2,3 
as a function of the frequency for a Cen A & B (models 
Abv & Bbv)- The "period" V of the small oscillation is 
203 //Hz for Abv and 255 ^Hz for Bbv 
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with frequency for radial orders n larger than about 16 for 
aCenA and 10 for aCenB. 

In the high frequency range the large and small fre- 
quency spacings which characterize the p-mode spectrum 
are usually estimated by analytical fits of the frequencies 
and of the small frequency separations. The numerical fre- 
quencies are fitted by the following polynomial expression 
(Berthomieu et al. 1993b| ): 



t I 
v n ,e = vat + Av e (n + - - n ) + a e (n + - - n ) 2 . 

The quantity 8i/£ t t+2 varies almost linearly with the fre- 
quency or the radial order. The small spacing is estimated 
by using the fit: 



5v L 



5vtj+2 + St(n - n ). 



As in the solar case, we consider a set of 9 modes cen- 
tered at no ~ 21, which corresponds approximately to the 
middle of the range of the expected excited frequencies 
(i.e. 1 — 3mHz) for a CenA, according to Houdek ( 199(f ). 
The fitted coefficients ae are very small. The quantity Ave 
does not much depend on the degree, so that Ave ~ Av 
and it characterizes the large frequency spacing. The small 
spacings are conventionally measured by 5vt,i+2- Table @ 
shows the quantities Vqq, Avq, Svq.2 and 8v\^ which have 
been computed for the models of a Cen A & B given in 
Table § 

These large and small separations of frequency, which 
characterize the p-mode oscillation spectrum, depend on 
the stellar mass and age and slightly decrease with the age 



(Christensen-Dalsgaard 1984, Audard et al. 1995). vqq and 



Ai^o are mainly related to the envelope structure of the 
stellar model and they vary proportionally to M^ii", 
while the values of 6vq,2 and <5^i ; 3 reflect the structure in 
the core. 

The computed values of Ai^o are respectively smaller 
(respt. larger) than the solar ones due to larger (respt. 
smaller) masses of a Cen A (respt. aCenB). These values 
do not depend much on the description of the convection 
or on the convective core overshoot, as can be seen from 
a comparison of the models Abv, A-CMi A ov and Bbv, 
Bcm, B ov . On the contrary, the value of Avq obtained 
for the oldest model Aqd of a CenA, which has also a 
larger mass, is significantly smaller by about 5 /xHz. The 
result is opposite for aCenB. All these differences can be 
accounted for by the differences in mass and radius. 

The small separations Svq,2 and Sv\^ decrease with 



stellar age and mass (Christensen-Dalsgaard 1984). The 
values of Svq^ and 5v\ ^ obtained for our calibrated mod- 
els are comparable to those of Pourbaix et al. (1999). As 



expected, they decrease with the age for the models with- 
out core overshoot, and thus are significantly smaller for 
the Aqd model of a Cen A. Future asteroseismic observa- 
tions could help to discriminate between such models. 

The lower panel of Table ^ presents the variations 
A[A^ ], A\8vq^\ and A\8vx .3] for the models computed 



with the MLTbv theory and extreme modeling parame- 
ters a, t a ccn, Y. and [§]i within the confidence domains 
presented in Table [l| The variations of Avq with a, i Q ccn 
and Yi are essentially due to the difference of model radii. 
The small spacings Svq^ and Svi.3 depend mainly on the 
age and are less sensitive to the differences in radius. 



According to Gough (1991), further information on the 



stellar structure can be provided by low degree oscilla- 
tions, from the second order difference of frequencies: 

A nt e = v nt i — 2v n+ i t ( + V n+ 21- 

When plotted as a function of the frequency (Figure ||), 
this quantity has a sinusoidal behavior with two differ- 
ent "periods" of order 210 //Hz and 700 /^Hz. The larger 
period has the largest amplitude and is due to the rapid 
variation of the adiabatic index Ti in the Hell ionization 
zone. Its contribution to the frequency is also clearly vis- 
ible in Fig. 0. The lower period V is an indication for a 
discontinuity in the derivative of the sound velocity at the 
base of the convection zone and is the inverse of twice the 
travel time of the sound from the surface to the base of 



the convection zone (e.g. Audard & Provost 1994) 



v- 1 



dr 



flea 



The predicted value of V is larger for a Cen B than for 
a Cen A, mainly due to a deeper convection zone for 
aCenB (see Table ||). The corresponding amplitude is 
sensitive to different processes, like convective penetration 



below the convection zone (e.g. Berthomieu et al. 1993a) 



8.3.2. g-mode oscillations 

In the low frequency range, the frequencies are mainly de- 
termined by the Brunt- Vaissala frequency N . The period 
of low degree gravity modes is proportional to a charac- 
teristic period: 

1 /" Rzc N 
v n ,£ ~ P (n + Eg) = — / — dr. 

ilT •/Rod r 

The integral, which defines Pq , is taken in the inner radia- 
tive zone, i.e. from the radius of the convective core R co for 
a Cen A, or from the center for a Cen B, to the base of the 
external convection zone R zc . For stars without convective 
core, like the Sun or a CenB, ee depends on the degree of 
the oscillation, eg = £/2 + e; for stars with convective core 
(like a Cen A) ee = e (Christensen-Dalsgaard 1984). Ta- 
ble gives the values of Pq for the models of a CenA & B. 
When the calibration of a, Cen A is made with a core over- 
shoot (which increases the extent of the convective core) 
Pq is significantly larger (by 23%). The other differences 
between the values of Pq for the models of a Cen A & B 
are accounted for by the differences in mass and in radius. 
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9. Discussion and conclusion 

Detailed evolutionary calculations of the visual binary 
a Centauri, including the pre main-sequence have been 
performed using the recent mass determinations of Pour- 
baix et al. ( 1999| ) . Models have been constructed using the 
CEFF equation of state, OPAL opacities, NACRE ther- 
monuclear reaction rates and microscopic diffusion. We 
have revisited the effective temperatures, surface gravi- 
ties and metallicities, using published spectroscopic data 
and taken these quantities as observational constraints. 
We have determined the most reliable solution within the 
confidence domains of the observable constraints via a 



X 2 minimization. Each solution is characterized by 



P 



\tot Cen ) Yi,[ a a, o-b}, where i Q Ccn is the age of the sys- 
tem, Y{ the initial helium content, [^]; the initial metal- 
licity and a\ and ae the convection parameters of each 
star model. We obtained calibrations using different con- 
vection theories and adapted values for the mixing-length 
parameter of each component. 

With the basic Bohm-Vitense ( |1958| ) mixing- 
length theory of convection, we derived pev = 
{2710 Myr, 0.284, 0.257, 1.53, 1.57}. With a con- 
vective core overshoot of 0.20 H p , we obtained 
p ov = {3530 Myr, 0.279 , 0.264 , 1.64, 1.66}. With the 
Canuto & Mazitclli (1991, 1992) convection theory, we 
get pom = {4086 Myr, 0.271, 0.264, 0.964, 0.986}. Using 
the mass values and the observational constrai nts o f 
the "best" model of Guenther & Demarque ( 2000| ), 
with the basic mixing-length theory, we obtained 
Pgd = {5638 Myr, 0.300, 0.296, 1.86, 1.97}. 

We have also performed \ 2 minimizations forcing the 
use of a unique value for the convection parameters of 
models of a Cen A & B. We have not reported the results 
as they are not significantly different from the former ones. 

The most striking fact in our results is the small val- 
ues obtained for the ages which are noticeably smaller 
than prev ious studies, except th ose of Boesgaard & Ha- 
gen (1974) and Pourbaix et al. (1999). With respect to 
the recent models of Guenther & Demarque ( [2000 ) this is 
due mainly to the mass discrepancies resulting from the 
small differences in distances. The ~ 5.5% increase be- 
tween the masses of Pourbaix et al. ( 1999|) we used and 
those used by Guenther & Demarque ( 2000| ) results from 
by the 1.5% disparity between the revised Hipparcos par- 
allax (Soderhjelm [2000 ) used by Guenther & Demarque 
and the orbital parallax determination of Pourbaix et al. 
As revealed by the x 2 minimization, in the case of the 
present calibration of a Cen, the determination of the age 
appears to be more sensitive to the mass differences than 
to the basic observational atmospheric constraints namely, 
effective temperature, gravity - or luminosity - and metal- 
licity. According to our results, no really satisfactory cri- 
terion allows to discriminate among different sets of mod- 
eling parameters: they all generate models which fit the 
observational constraints based on H-R diagram analysis 



and metallicities. A naive argument, however, can plead 
in favor of a subsolar age: as the a Cen binary system is 
formed of metal enriched material, it formed from more 
processed interstellar matter than the Sun; this may indi- 
cate that a Cen is younger than the solar system. 

For the models computed with the basic MLTbv 
theory and fitted on our observational constraints, the 
mixing-length parameters are noticeably smaller than the 
convection parameter a© of the solar model calibrated 
with the same physics; the diff erence s are larger than ex- 
pected from the Ludwig et al. ( 1999|) simulations. For the 
models fitted to the observational const raints of the "best" 
model of Guenther & Demarque ( [2000 ), the values derived 
for a bracket the calibrated solar value. For all models 
computed with MLTbv theory we note that the smaller 
the age, the smaller the mixing-length parameters. 

For models computed with MLTcm convection theory, 
in accordance with Canuto & Mazitelli (1991, 1992) we 



obtain mixing length parameters both close to unity and 
to the convection parameter a© of the solar model cali- 
brated with the same physics, a a — Qb — a© — L This 
indicates that the convection theory of Canuto & Mazitclli 
may support the assumption of a universal convection pa- 
rameter and it seems to provide better fits to observations, 
as already pointed out by helioseismology as aforesaid in 
Sect.|. 

For the a Cen A models fitted on our observational 
constraints, as in the solar model, microscopic diffusion 
alone is not efficient enough to account for the observed 
surface lithium depletion. This may indicate that an un- 
known physical process, at work beneath the outer convec- 
tion zone, reinforces the microscopic diffusion and gravi- 
tational settling to transport the material down to the 
lithium burning zone. For aCenB, the only available ob- 
servation of the surface lithium abundance is an upper 
limit (Chmielewski et al. 1992): all aCenB models pre- 



dict values compatible with this observation except model 
Bbv- The surface lithium depletion of models calibrated 
with the Guenther & Demarque ( 2000] ) observational con- 
straints are close to the observations and, in that sense, 
they appear to be more satisfactory than models fitted 
on our observational constraints: in these models there is 
no need for additional physical process to transport the 
lithium to its burning zone. But this imply that the small 
lithium depletion found in the standard solar model and 
the too large efficiency of the microscopic diffusion in more 
massive stars result from physical processes which are not 
active in a Cen, though masses and rotation status are 
close to solar ones. 

All our models with large masses have initial helium 
contents close to the solar value Y© — 0.003 < Yi < Yq + 
0.010, while the models with low masses have a high initial 
helium content Yi = Yq +0.026. With a primordial helium 
abundance of Yq « 0.235 we get a galactic enrichment of 
« 2.0, both for the Sun and low mass models and 
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^2 = 1-2 — 1-6 for high mass models. The differences in 
^2 between our a Cen models and the Sun are compatible 
with the scattering found in the solar neighborhood (Pagel 
& Portinari 1998) and in other binary system calibrations 
(Fernandes et al. 1998). 

We have computed the large and small frequency spac- 
ings of acoustic oscillations for all the models. The large 
separation Avq for the three models with large masses are 
within 2 /iHz, despite their differences in age and physics. 
They are such that Ai/q ~ 107 — 108 /iHz. The variations 
of A^o taking into account the uncertainties in the ob- 
servable constraints, effective temperatures, gravities and 
metallicities, are of about ±1 /iHz. For the model Aqd 
of small mass, the large separation, Aisq ~ 102 /iHz, is in 



agreement with the value of Guenther & Demarque ( 2000) 



and significantly lower for our other models. The differ- 
ences in Avq are mainly explained by the differences in 
mass and radius. 

For the three models with large masses the small sepa- 
rations (5^o,2 vary from 7.5 to 9.1 /iHz for a Cen A and are 
of about 12 /iHz for a Cen B. For each model the variations 
in (5^o,2 within the observable uncertainties are ±1 /iHz. 
For the model Aqd, <5^o,2 is much smaller and in agree- 
ment with the value of Guenther & Demarque. All these 
differences in 5i/o,2 are mainly related to the differences in 
central hydrogen content, hence in the age. 

These results show that the determination of Avq and 
$vo,2 by seismological observations would help to discrim- 
inate between the models of a Cen A computed with dif- 
ferent masses and to confirm or not the new determination 
of the masses by Pourbaix et al. This implies an improve- 
ment of the accuracy of the observables used to constrain 
the calibrated models of a Cen A & B. Concerning the 
comparison with the seismic observations, the large split- 
ting estimated by Pottasch et al. (|1992| ) and Edmonds & 
Cram fll995| ) favor Pourbaix et al fll999| ). masses The dif- 
ferent possible estimations for the large and small spacings 
Avq and <5z/o,2 by Kjeldsen et al. (1999a) do not allow to 
discriminate between the models. We note, however, that 
their estimation Auq = 100.8 /iHz and <5^o.2 = H-7 /iHz is 
highly improbable. 

We conclude that, even for a Cen, the best known bi- 
nary system, the models are not strongly enough con- 
strained by the available astrometric, photometric and 
spectroscopic data. In order to deeply test stellar physics 
additional information on the internal structure is needed. 
Up to now ground-based observations give tentative evi- 
dence for acoustic oscillations in aCen A. In a few years 
from now, one can expect that asteroseismology from 
space e.g. COROT (Baglin et al. ggfgb, MONS (Kjeld- 
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Table 5. All symbols have their ordinary meaning. {The brackets indicate values derived from basic formula.} 
Top, left panel: adopted observational data to be reached by the calibration. 



Top, right panel: Observational constraints of the "best" model of Guenther & Demarque (200C). 
Bottom panel: Characteristics of a Cen A & B and Sun models computed with the same input physics. The aCenA 
(respt. aCenB) models are named "A.." (respt. "B.."). The first four rows recall some items of Table |l|: the ages 
(Myr), the initial helium mass fractions, the initial heavy element to hydrogen mass ratios, and the mixing- length 
parameters. The next three rows give the effective temperatures in K, the surface gravities and the metallicitics. 
The six next rows present respectively the luminosities, the total radii in solar units, the surface mass fractions of 
hydrogen and helium, the ratios of heavy element to hydrogen and the lithium abundances. In the next rows R cz and 
T cz are respectively the radius and the temperature (in MK) at the base of the external convection zone, R co is the 
radius of the convective core (including overshoot for model A ov & B ov ). At center, T c , p c , X c , Y c are respectively the 
temperature (in MK), the density (in gem -3 ), the hydrogen and the helium mass fractions. 



Models bv, ov & CM 
a Cen A a Cen B 



Models GD 
a Cen A a Cen B 



M/M & 

T e ff 

logs 

L/L e 



1.16 ±0.031 0.97 ±0.032 1.1015 ± 0.008 0.9159 ±0.007 

5790 ± 30 K 5260 ± 50 K 5770 ± 50 K 5300 ± 50 K 

4.32 ± 0.05 4.51 ±0.08 {4.28 ±0.02 4.54 ±0.03} 

0.20 ±0.02 0.23 ±0.03 0.22 ± 0.02 0.26 ± 0.04 

{1.534 ±0.103 0.564 ±0.064} 1.572 ±0.135 0.509 ± 0.06 



Models with MLTbv theory 



Models with MLT C m theory 



A BV 


Bbv 




B ov 


A GD 


Bgd 


0BV 


Aqm 


Bcm 


©CM 


2710 


3530 


5640 


4685 


4086 


4685 


0.284 


0.279 


0.300 


0.274 


0.271 


0.274 


0.0443 


0.0450 


0.0480 


0.0279 


0.0450 


0.0279 


1.53 


1.57 


1.64 


1.66 


1.86 


1.97 


1.93 


0.96 


0.99 


1.03 


5795 


5269 


5795 


5269 


5759 


5332 


5778 


5794 


5268 


5778 


4.324 


4.508 


4.316 


4.509 


4.286 


4.523 


4.438 


4.324 


4.522 


4.438 


0.20 


0.23 


0.20 


0.23 


0.22 


0.25 


0.000 


0.20 


0.23 


0.000 


1.527 


0.571 


1.555 


0.569 


1.565 


0.547 


1.000 


1.526 


0.552 


1.000 


1.228 


0.909 


1.239 


0.907 


1.259 


0.867 


1.000 


1.228 


0.893 


1.000 


0.722 


0.701 


0.729 


0.708 


0.711 


0.695 


0.737 


0.735 


0.717 


0.737 
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0.270 


0.243 


0.262 
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0.275 


0.245 


0.236 


0.253 


0.245 


2.72 


1.10 


2.54 
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1.94 


-3.32 


2.46 


2.33 


-0.54 


2.46 
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0.0385 
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0.0433 


0.0245 
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0.0245 
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0.698 


0.758 
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0.707 


0.675 


0.715 


0.741 


0.685 


0.715 
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2.585 


1.644 


2.653 


2.020 


2.846 


2.178 


1.801 


2.761 


2.178 


0.026 




0.054 




0.061 






0.038 






17.00 


13.66 


17.46 


13.76 


19.43 


14.26 


15.77 


17.74 


13.71 


15.76 


126.0 


96.28 


125.5 


100.8 


172.0 


123.7 


155.9 


146.8 


102.4 


155.9 


0.374 


0.540 


0.387 


0.510 


0.178 


0.400 


0.330 


0.283 


0.497 


0.330 


0.594 


0.428 


0.580 


0.458 


0.788 


0.566 


0.649 


0.684 


0.470 


0.649 
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(#)• 
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Table 6. Partial derivatives of observables logT e ffA, log 7^-, [tt]a' l°S^effB, «b, log and [§] B with respect to 
modeling parameters i a Ccn(Myr), a a, Yi and [^] of models Abv & Bbv- 



aiogT offA Slog^ 9[f] A aiogT offB 91og^| £>[§]. 



9ic« Con 


-2.15 


4.16 10~ 5 


-2.30 10~ 5 


1.63 10~ 6 


2.32 10~ 5 


-9.2410" 


<9o?a 


0.0380 


8.23 10~ 3 


0.0541 


0. 


0. 


0. 


da B 


0. 


0. 


0. 


0.0306 


0.0213 


5.67 10" 3 


dYi 


0.469 


3.85 


-0.911 


0.621 


3.55 


-0.185 




-2.12 10 -3 


-0.014 


0.024 


-2.40 10 -3 


-0.0129 


0.0205 
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Table 8. Low degree frequencies for aCenA and B. 







aCenA, model Abv 






aCenB, model Bbv 




n 


£ = 


1=1 


I = 2 


I = 3 


£ = 


I = 1 


£ = 2 


£ = 3 


1 


216.50 


239.04 


332.72 


354.66 


280.33 


313.25 


404.94 


460.11 


2 


330.58 


357.01 


414.45 


450.61 


439.72 


504.26 


575.00 


635.19 


3 


433.74 


473.41 


525.29 


567.19 


594.94 


674.45 


748.46 


810.92 


4 


539.84 


589.81 


640.82 


683.00 


765.25 


845.17 


916.83 


980.26 


5 


651.62 


703.50 


754.02 


797.53 


932.81 


1013.63 


1085.89 


1148.42 


6 


763.65 


817.42 


867.89 


911.05 


1101.83 


1179.48 


1251.72 


1315.75 


7 


876.96 


929.96 


980.65 


1024.62 


1268.49 


1346.02 


1416.99 


1480.67 


8 


989.90 


1043.11 


1093.39 


1136.78 


1433.70 


1510.26 


1581.42 


1645.03 


9 


1102.65 


1154.54 


1204.33 


1247.43 


1598.27 


1673.35 


1742.83 


1806.05 


10 


1213.79 


1264.47 


1312.49 


1354.49 


1759.39 


1833.61 


1901.74 


1963.56 


11 


1321.79 


1370.95 


1418.40 


1460.41 


1917.69 


1989.93 


2057.43 


2119.37 


12 


1427.80 


1476.58 


1523.87 


1566.71 


2073.10 


2145.00 


2212.06 


2274.53 


13 


1533.43 


1583.08 


1631.25 


1674.60 


2227.40 


2299.89 


2367.70 


2430.43 


14 


1640.95 


1690.78 


1738.77 


1782.29 


2382.75 


2455.04 


2522.86 


2585.94 


15 


1748.45 


1798.25 


1845.79 


1888.78 


2537.58 


2609.85 


2677.04 


2739.71 


16 


1855.24 


1904.34 


1951.62 


1994.97 


2691.18 


2762.91 


2830.18 


2892.90 


17 


1961.01 


2010.40 


2057.91 


2101.70 


2843.92 


2915.35 


2982.68 


3046.16 


18 


2067.20 


2117.07 


2165.32 


2209.86 


2996.13 


3068.24 


3136.07 


3200.00 


19 


2174.62 


2225.10 


2273.55 


2318.42 


3149.16 


3221.60 


3290.20 


3354.77 


20 


2282.71 


2333.49 


2382.26 


2427.33 


3303.04 


3375.80 


3444.50 


3509.53 


21 


2391.27 


2442.08 


2490.76 


2536.02 


3456.99 


3530.11 


3599.13 


3664.29 


22 


2499.56 


2550.55 


2599.41 


2644.81 


3611.26 


3684.29 


3753.67 


3819.43 


23 


2608.02 


2659.04 


2708.02 


2753.86 


3765.55 


3838.97 


3908.55 


3974.76 


24 


2716.48 


2767.90 


2817.19 


2863.25 


3920.10 


3993.88 


4064.07 


4130.66 


25 


2825.46 


2877.03 


2926.58 


2973.04 


4075.36 


4149.28 


4219.77 


4286.93 


26 


2934.69 


2986.58 


3036.24 


3082.83 


4230.79 


4305.15 


4375.88 


4443.28 


27 


3044.15 


3096.07 


3145.94 


3192.80 


4386.60 


4461.07 


4532.20 


4600.01 


28 


3153.68 


3205.77 


3255.68 


3302.66 


4542.70 


4617.33 


4688.59 


4756.82 


29 


3263.18 


3315.36 


3365.50 


3412.70 


4698.84 


4773.76 


4845.36 


4913.76 


30 


3372.84 


3425.10 


3475.26 


3522.63 


4855.34 


4930.30 


5002.17 


5070.99 


31 


3482.42 


3534.76 


3585.10 


3632.58 


5011.94 


5087.14 


5159.10 


5228.14 


32 


3592.07 


3644.40 


3694.70 


3742.31 


5168.62 


5243.95 


5316.22 


5385.43 


33 










5325.52 


5400.85 


5473.22 


5542.77 



